The nematic phase of a system of long hard rods* 
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We consider a two-dimensional lattice model for liquid crystals consisting of long 
rods interacting via purely hard core interactions, with two allowed orientations 
defined by the underlying lattice. We rigorously prove the existence of a nematic 
phase, i.e., we show that at intermediate densities the system exhibits orientational 
order, either horizontal or vertical, but no positional order. The proof is based on a 
two-scales cluster expansion: we first coarse grain the system on a scale comparable 
with the rods' length; then we express the resulting effective theory as a contour's 
model, which can be treated by Pirogov-Sinai methods. 

Dedicated to the 70th birthday of Giovanni Gallavotti 

1. INTRODUCTION 

In 1949, L. Onsager [30] proposed a statistical theory for a system of elon- 
gated molecules interacting via repulsive short-range forces, based on an explicit 
computation of the first few Mayer's coefficients for the pressure. Onsager's the- 
ory predicted the existence at intermediate densities of a nematic liquid crystal 
phase, that is a phase in which the distribution of orientations of the particles is 
anisotropic, while the distribution of the particles in space is homogeneous and 
does not exhibit the periodic variation of densities that characterizes solid crys- 
tals (periodicity in all space dimensions) or smectic liquid crystals (periodicity in 
one dimension). 

From a microscopic point of view, the most natural lattice model describing 
elongated molecules with short-range repulsive forces is a system of rods of length 
k and thickness 1 at fixed density p (here p = average number of rods per unit 
volume), arranged on a cubic lattice, say a large squared box portion of Z^, 
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and interacting via a purely hard core potential. Even though very natural, this 
model is not easy to treat and its phase diagram in the plane (p, k) is still not 
understood in many physically relevant parameters' ranges. Of course, for all 
fc's, at very small density there is a unique isotropic Gibbs state, invariant under 
translations and under discrete rotations of 90°; this can be proved by standard 
cluster expansion methods. If A; = 2, it is known [17] that the state is analytic 
and, therefore, there is no phase transition, for all densities but, possibly, at the 
close packing density, i.e., at the maximal possible density Pmax = ^/k. If /c is 
sufficiently large {k > 7 should be enough [13]) there is numerical evidence [13, 26] 
for two phase transitions as p is increased from zero to the maximal density. The 
first, isotropic to nematic, seems to take place at a p'^c^ ~ Ci/k'^, while the second, 
nematic to isotropic, seems to take place at pi^^ ~ p^nax — C2/k^. These findings 
renovated the interest of the condensed matter community in the phase diagram 
of long hard rod systems and stimulated more systematic numerical studies of 
the nature of the critical points at p^c^ and pi^^ [8, 9, 24, 25, 27, 28]. Prom a 
mathematical point of view there is no rigorous proof of any of these behaviors 
yet, with the exception of the "trivial" case of very low densities: namely, there 
is neither a proof of nematic order at intermediate densities, nor a proof of the 
absence of orientational order at very high densities, nor a rigorous understanding 
of the nature of the transitions. 

In this work wc give a rigorous proof to some of the conjectures stated above 
on the nature of the phase diagram of long hard rods systems. More precisely, 
we show that well inside the interval {p^c\ Pc^), the system is in a nematic phase, 
i.e., in a phase characterized by two distinguished Gibbs states, with different 
orientational order, either horizontal or vertical, but with no positional order. To 
the best of our knowledge, this is the first proof of the existence of a nematic phase 
in a microscopic model with molecules of fixed finite length and finite thickness, 
interacting via a purely repulsive potential. In this respect, our result is a strong 
confirmation of Onsager's proposal that orientational ordering can be explained 
as an excluded volume effect. 

Our proof is based on a two-scales cluster expansion method, in which wc first 
coarse grain the system on scale k; we next realize that the resulting effective 
model can be expressed as a contour model, reminiscent of the contour theory 
for the Ising model at low temperatures. However, contrary to the Ising case, the 
contour theory we have to deal with here is not invariant under a 7? symmetry: 
therefore, wc cannot apply the Peierls' argument and we need to make use of a 
Pirogov-Sinai method. 

Of course, our proof leaves many questions about the phase diagram of long 
hard rod systems open, the most urgent being, we believe, the question about 
the nature of the densely packed phase at p > pi^^: can one prove the absence 
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of orientational order, at least at close packing? Is the densely packed phase 
characterized by some "hidden" (striped-like) order? Progress on these problems 
would be important for the understanding of the emergence of hidden order in 
more complicated systems than elongated molecules with purely hard core inter- 
actions, in which short range repulsion competes with attractive forces acting on 
much longer length scales. 

Previous results. There is a limited number of papers where important previous 
results on the existence of orientational order in lattice or continuum models for 
liquid crystals were obtained, related to the ones found in this work. 

A first class of liquid crystal models that has been considered in the literature 
describe long rods with purely repulsive interactions and discrete orientations, like 
ours; of course, the case of continuous orientations would be of great interest, but 
its treatment appears to be beyond the current state of the art. In [16, 19, 22], the 
existence of orientational order for different variants of lattice gases of anisotropic 
molecules with repulsive interactions was proved, by using Peierls-like estimates 
and cluster expansion; however, in all these cases, orientational comes together 
with translational order, which is not the case in a nematic phase. A continuum 
version of the model in [22], i.e., a continuum system of infinitely thin rods 
with two allowed orientations and hard core interactions, was later proved to 
have a phase transition from an isotropic to a nematic phase [5, 33], by using 
improved estimates on the contours' probabilities and a Pirogov-Sinai method. 
More recently, the existence of an isotropic to nematic transition in an integrable 
model of polydisperse long rods in 7? with hard core interactions was proved [20] , 
by mapping the partition function of the polydisperse hard rods gas into that of 
the nearest neighbor 2D Ising model. 

A second class of liquid crystal models studied in the literature assumes the 
existence of attractive forces favoring the alignment of the molecular axes: in fact, 
in some cases, the attraction is expected to originate from the inter-molecular 
Coulomb interaction [29] and to play a more prominent role than the Onsager's 
excluded volume effect. The emergence of a nematic phase in such models was 
first understood at the mean field level [3, 7, 29, 31]. Later, it was understood 
that in the presence of attractive forces, even the monomer-dimer system can 
exhibit an oriented phase at low temperatures, as proved in [18] by refiection 
positivity methods; the absence of positional order for the same model, known as 
the Heilmann-Lieb's model [18], was then proved on the basis of cluster expan- 
sion methods [23]. Remarkably, if attractive forces favoring the alignment of the 
molecules' axes are allowed, there are models displaying a full 0{m) orientational 
symmetry, m > 2, for which it is possible to rigorously prove the existence of 
nematic order (or quasi- long range order, depending on the dimensionality). In 
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particular, in [1, 2, 35] certain d-dimensional lattice-gas models describing par- 
ticles with an internal ("spin") continuous orientational degree of freedom were 
introduced; the existence of orientational order was proved both in > 2 for 
short-ranged interactions and in = 1,2 with sufficiently long-ranged interac- 
tions, via a combination of infrared bounds and chessboard estimates [10]. In [15], 
a proof of the existence of orientational quasi long range order a'la Kosterlitz- 
Thouless was given for a similar system in = 2 with short ranged interactions, 
by using a combination of the Gruber-Griffiths method [14], originally applied to 
the study of an orientational phase transition in a continuum system of particles 
with internal Ising-hke degrees of freedom, and of the Prohlich-Spencer method 
[11], originally apphed to the study of the Kosterhtz-Thouless transition in the 
classical two-dimensional XY model. 

Summary. The rest of the paper is organized as follows. In Section 2 we 
"informally" introduce the model, state the main results and explain the key 
ideas involved in the proof. In Section 3 we define the model and state the main 
theorem (Theorem 1 below) in a mathematically precise form. In the following 
sections we prove Theorem 1: in Section 4 we rewrite the partition function with 
q boundary conditions in terms of a sum over contours' configurations, where 
the contours are defined in a way suitable for later application of a Pirogov-Sinai 
argument. In Section 6 we prove the convergence of the cluster expansion for 
the pressure, under the assumption that the activity of the contours is small and 
decays sufficiently fast in the contour's size. In Section 7 we complete the proof 
of convergence of the cluster expansion for the pressure, by inductively proving 
the desired bound on the activity of the contours. Finally, in Section 8 we adapt 
our expansion to the computation of correlation functions and we prove Theorem 
1. 



2. THE MODEL 

We consider a finite square box A C of side L, to be eventually sent to 
infinity. We fix k and the average density p G (0, 1/A;). The finite volume Gibbs 
measure at activity z gives weight 2;" to every allowed configuration of n rods: 
we say that a configuration is allowed if no pair of rods overlaps. Of course, one 
also needs to specify boundary conditions: we consider, say, periodic boundary 
conditions, open boundary conditions, horizontal or vertical boundary conditions, 
the latter meaning that all the rods within a distance ~ k from the boundary of 
A are horizontal or vertical - see below for a more precise definition. The grand 
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canonical partition function is: 

Za(^)=^;^«^^ (2.1) 



n>0 



where is the number of allowed configurations of n rods in the box A, in 
the presence of the prescribed boundary conditions. Note that = for all 
n > \A\/k, which shows that Z\{z) is a finite (and, therefore, well defined) sum 
for all finite A's. The activity z is fixed in such a way that 
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lim ^= lim ^^^lifi— ^ = (2.2) 

|A|->oo |A| |A|-^oo |A| Zj^{z) 

The goal is to understand the properties of the partition function and of the 
associated Gibbs state in the hmit |A| — > oo at fixed p. An informal statement 
of our main result is the following. 

Main result. For k large enough, if k^^ ^ p ^ k"^ , the system admits two 
distinct infinite volume Gibbs states, characterized by long range orientational 
order (either horizontal or vertical) and no translational order, selected by the 
boundary conditions. 

Sketch of the proof. The idea is to coarse grain A in squares of side i ^ k/2. 
Each square is large, since in average it contains many (~ pk"^ ^ 1) rods. On 
the other hand, its side C. is small enough to ensure that only rods of the same 
orientation are allowed to have centers in the same square. This means that the 
partition function restricted to a single square contains only sums over vertical 
or horizontal configurations. Let us consider the case where the rods are all 
horizontal (vertical is treated in the same way) . A typical horizontal configuration 
consists of many (~ pk"^) horizontal rods with centers distributed approximately 
uniformly (Poisson-like) in the square, since their interaction, once we prescribe 
their direction, is very weak: they "just" have a hard core repulsion that prevents 
two rods to occupy the same row, an event that is very rare, since the density 
of occupied rows (~ pk'^/k) is very small, thanks to the condition that p -C l/k. 
Because of this small density of occupied rows, we are able to quantify via cluster 
expansion methods how close to Poissonian is the distribution of the centers in 
the given square (once we condition with respect to a prescribed orientation of 
the rods). 

To control the interaction between different squares we use a Pirogov-Sinai 
argument. Each square can be of three types: (i) either it is of type +1, if it 
contains only horizontal rods, (ii) or it is of type —1, if it contains only vertical 
rods, (iii) or it is of type 0, if it is empty. The values —1, 0, +1 associated to each 
square play the role of spin values associated to the coarse grained system. The 
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interaction between the spins is only finite range and squares with vertical (+1) 
and horizontal (—1) spin have a strong repulsive interaction, due to the hard core 
constraint. On the other hand, the vacuum configurations (the spins equal to 0) 
are very unlikely, since the probability of having a large deviation event such that 
a square of side I is empty is expected to be exponentially small ~ exp{— cp/c^}, 
for a suitable constant c. 

Therefore the typical spin configurations consist of big connected clusters of 
"uniformly magnetized spins", cither of type +1 or of type —1 separated by 
boundary layers (the contours), which contain zeros or pairs of neighboring op- 
posite spins. These contours can be shown to satisfy a Peierls' condition, i.e., 
the probability that a given contour occurs is exponentially small in the size 
of its geometric support. The contour theory is not symmetric under spin fiip 
and, therefore, we are forced to study it by the (non-trivial although standard) 
methods first introduced by Pirogov and Sinai [32] . 

Before we move to discuss the details of our proof, let us state our main results 
in a mathematically more sound form. 



Definitions. For any region X C we call Vtx the set of rod configurations 
R = {ri, . . . ,r„} where all the rods belong to the region X. A rod r "belongs 
to" a region X if the center of the rod is inside the region, in which case we 
write r E X. Here each rod is identified with a sequence of k adjacent sites 
of in the horizontal or vertical direction. If k is odd, the center of the rod 
belongs to the lattice Z^ itself and, therefore, the notion of "rod belonging to X" 
is unambiguously defined. On the contrary, if k is even, the geometrical center 
of the rod does not belong to the original lattice Z^; however, for what follows, 
it is convenient to pick one of the sites belonging to r and elect it to the role of 
"center of the rod" : if r is horizontal (vertical) , we decide that the "center of r" 
is the site of r that is closest to its geometrical center from the left (bottom). 
We shall also say that: a rod r ''touches" a region X, if r fl X 7^ 0; a rod r "is 
contained in" a region X, if r fl X'^ = 0, in which case we write r C. X. 

The rod configurations in Qx can contain overlapping and even coinciding 
rods; we denote by R{r) the multiplicity of r in i? G Qx- The grand canonical 
partition function in X with open boundary conditions is 



3. 



MAIN RESULTS 




(3.1) 



Renx 
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where \R\ :— R{r) and (p{R) implements the hard core interaction: 

V.(fl)= n V(r,r'), ^M = \lZn''r-Zl. ''''^ 

r,r'GR y ' 

Let i := [/c/2'1 and assume that A C is a square box of side divisible by 4£. 
We pave A by squares of side called "ti/es", and by squares of side 4£, called 
'''' smoothing squares'^ The lattice of the tiles' centers is a coarse grained lattice 
of mesh £, called A'; similarly, the lattice of the smoothing squares' centers is 
a coarse grained lattice of mesh 4£, called A". Given ^ G A', the tile centered 
at ^ is denoted by A^; given a e A", the smoothing square centered at a is 
denoted by Sa- Given two sets X,Y C A, we indicate their euchdean distance 
by dist{X,Y) — mmxex,yeY\x — y\- li X and Y are union of tiles, we shall 
also indicate by X',Y' C A' the coarse versions of X and Y, i.e., the sets of 
sites in A' such that X — U^gx'A^ and Y — U^gyA^. The distance between X' 
and Y' is denoted by dist{X' ,Y') and their rescaled distance by dist'{X' ,Y') :— 
£~^dist{X' ,Y'); with these conventions, if and 77 are nearest neighbor sites on 
A', then dist(^, 77) = I^C ~ ^| = ^ and dist'(^,r7) = 1. The complement of A is 
denoted by Ac := \ A and its coarse version by A^, with obvious meaning. 

The size of the tiles is small enough to ensure that if one vertical (horizontal) 
rod belongs to a given tile, then all other rods belonging to the same tile and 
respecting the hard core repulsion condition must be vertical (horizontal). If a 
tile is empty, i.e., no rod belongs to it, then we assign it an extra fictitious label, 
which can take three possible values, either or + or — . A rod configuration 
R & Q\ (combined with an assignment of these extra fictitious labels) induces a 
spin configuration a = {cr^j^gA' on A', G {—1,0, +1}, via the following rules: 

- (7^ = +1, if all rods belonging to are horizontal or if the tile is empty 
with the extra label equal to +, 

- (7^ = —1, if all rods belonging to A^ are vertical or if the tile is empty with 
the extra label equal to — , 

- cr^ 0, if A^ is empty with the extra label equal to 0. 

The corresponding set of rod configurations in the tile A^ is denoted by ^'2^'- 
(^Aj) configurations in Aj consisting either of horizontal 

(vertical) rods or of the empty configuration; similarly, fi^^ consists only of the 
empty configuration. 

Note that the grand canonical partition function in A with open boundary 
conditions can be rewritten as 

^o(A) = E E ' (3-3) 
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where ©a' := {— 1,0,+1}^' and 0,\{a) :— \J^^a'^'2^- Moreover, 

m--= [U<(0]¥>{R), (3.4) 
where the activity of a tile is defined as 

1—1 II = 0. 

The sign —1 is necessary to avoid over-counting of the empty configurations. Note 
that (f{R) depends both on a and on R; however, in order not to overwhelm the 
notation, we shall drop the label a. 

The partition function with q boundary conditions, q — ±, denoted by Z{A\q), 
can be defined in a similar fashion: 

m<i) = E E ^(^) (3-6) 

where 0^, C ©a' is the set of spin configurations such that dist'(,^, A'^) < 5 ^ 
= q. The number 5 appearing here is related to the choice of smoothing 
squares of side 4£: in fact, the condition that all the spins with dist'(,^, A'^) < 5 
are equal to q guarantees that all the smoothing squares adjacent to the boundary 
of A are uniformly "magnetized" with magnetization q and that, moreover, all 
such smoothing squares are surrounded by a 1-tile-thick peel of spins equal to 
q. These two conditions are convenient for an explicit construction of a contour 
representation for Z{A\q), as we will show below. 

Correspondingly, the ensemble (•)^ with q boundary conditions is defined by 

^^^^^^ZfL) ^ ^ ^(i^)^x(i?), (3.7) 

where Ax is a local observable, depending only on the restriction Rx of the rod 
configuration to a given finite subset X C A. The infinite volume states 
with q boundary conditions are defined by 

{Axr^lim {Ax)l, (3.8) 

|A|->-oo 

if the limit exists for all local observables Ax, X dl?. Our main results can be 
stated as follows. 

Theorem 1 If zk and {zk^)~^ are small enough, then the two infinite volume 
states {-y, q — ±, exist. They are translationally invariant and are different 
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among each other. In particular, ifx^o the projection onto the rod configurations 
such that e Q^^^ , then 

(X,7)' < ^-''"^ ' (3-9) 
for a suitable constant c. Moreover, let n^Q be the indicator function that is equal 
to 1 if a rod has a center in Xq G and otherwise, then 

P = (^..o)^ = Ko)" = ■^(l + 0{zk, e-^^'^')) (3.10) 

and 

p{x y) = {n^uy)^ = {n^uy}- =p'{l + 0(e-'^l--^l/^)) , (3.11) 
for a suitable c > 0. 

Eq.(3.9) proves the existence of orientational order in the system. Eqs.(3.10)- 
(3.11) prove the absence of translational symmetry breaking. These two behavior 
together prove that the system is in a nematic hquid crystal phase, as announced 
in the introduction. The rest of the paper is devoted to the proof of Theorem 
1, which is based on a two-scales cluster expansion. As it will be clear from the 
discussion in the next sections, our construction proves much more than what is 
explicitly stated in Theorem 1, namely it allows us to compute the averages of 
all the local observables in terms of an explicit exponentially convergent series. 

4. THE CONTOUR THEORY. 

The proof of the Theorem 1 will be spht in several steps. We start by develop- 
ing a representation of the partition function Z{A\q) with q boundary conditions 
in terms of a set of interacting contours. Later, we will adapt the contour expan- 
sion to the computation of the correlations. The contour theory can be studied by 
an adaptation of Pirogov-Sinai's method to the present context. See [32] for the 
original version of this method and [4, 21, 36, 37] for several alternative simpli- 
fied versions of it. In the following we will try to be as self-consistent as possible 
and to keep things simple, by avoiding as much as we can general and abstract 
settings. We first need some more definitions. 

Definition 1: sampling squares. Given a spin configuration a G 0^/, this 
induces a partition of A' into regions where the spins are "uniformly magnetized 
up or down" (i.e., regions where the spins are constantly equal to +1 or to —1) 
and boundary regions separating the "uniformly magnetized regions" among each 
other, which can possibly contain spins equal to zero. To make this more precise 
we introduce the notion of "sampling squares", defined as follows: given ^ G A', 
the sampling square associated to ^ is defined as = U^gA/; o<»7i-5i<^ ^r;, where 
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tile A 



sampling square S 




smoothing square S 



a 



FIG. 1. An example of tile, sampling square and smoothing square 



rji and r^j, i = 1,2, arc the coordinates of ^,?7 G A'. Note that if dist'(,^, A'^) > 1, 
then contains exactly 4 tiles. See Fig.l for an example. We say that a sampling 
square is 

• good if the spins inside are all equal cither to +1 or to —1. Each good 
sampling square comes with a magnetization m — ±1. 

• bad otherwise; note that each bad sampling square is such that either it 
contains at least one spin equal to zero, or it contains at least one pair of 
neighboring spins with opposite values, +1 and —1. 

Definition 2: connectedness, good and bad regions. Given a configu- 
ration a e 0A', we call 

5(a) = U^eA': (4.1) 

5{ is bad 

the union of all bad samphng squares. The "smoothening" of B{a) on scale M is 
defined as: 

B{a) = U aeA": Sa , (4.2) 

where the lattice A" and the smoothing squares Sa were defined in the paragraph 
following equation (3.2). 

Let X C A be a union of tiles: we say that X is connected if, given any pair 
of points x,y & X, there exists a sequence {xq = x,xi,. . . ,Xn-i,Xn = y) such 
that Xi & X and |xj — Xi_i\ = 1, for alH = 1, . . . , n. We also say that X is D- 
connected (with the prefix "D" meaning "diagonal") if, given any pair of points 
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x,y & X, there exists a sequence {xq — x,xi, . . . , Xn-i, Xn — y) such that Xi & X 
and \xi — Xi-i\ < \/2, for alH = 1, . . . , n (here |x — y| is the euchdean distance 
between x and y). 

The maximal D-connected components of B{(t) are denoted by Tj and arc the 
geometric supports of the contours that we will introduce below. The complement 
of the bad region, 

G{a) := A\B{a) , (4.3) 

can be split into uniformly magnetized disconnected regions, each of which is a 
union of tiles; these are denoted by Yj and ruj are the corresponding magnetiza- 
tions. 

Remcirks. 

1. Note that distinct D-disconnected bad regions in B{a), rj(o"), rj/((T) with 
j ^ j', do not interact directly; i.e., ip{R^,Rri) = 1 for all ^ G Tj, rj G Tj/. 
This is because I'jicr) and I'j'icr) are separated by at least one smoothing 
square (hence 4 tiles). Similarly, distinct uniformly magnetized discon- 
nected regions, Yj{a),Yji{a) G G{a) with j ^ f and magnetizations rrij, 
rrij/, do not interact directly; i.e., (fi{R^, Rq) — 1 for all ^ G Ij, 77 G Yji and 
for all i?^ G Q!^\ Rjj G Q,^^' . In fact, note that R^ and R^, can interact only 
in one of the following two cases: ^ and 77 are on the same row (column) 
and \C - vl < 2£> or |^i - 771I = |^2 - r]2\ = ^- If ^ e Yj and 77 G Yj> with 
j 7^ j', then the first case can occur only if |^ — 77] > 5i (in the horizontal 
or vertical directions, Yj and Yji are separated by at least one smoothing 
square), in which case R^ and R^ certainly do not interact, whatever is 
the alignment of the rods. In the second case necessarily rrij = rrij/, other- 
wise the sampling square containing both ^ and 77 would be bad and both 
tiles would belong to B{a) instead of G((t). Now, if rrij = rriji the rods in 
R^ have the same orientation as those in i?^, while their centers belong to 
different rows and columns and, therefore, do not interact. 

2. In terms of the definitions above, the set 0^^, C Ga' of spin configurations 
with q boundary conditions can be thought as the set of spin configurations 
such that all the contours' supports Tj C B{a) are D-disconnected from A'^ 
and separated from it by at least one smoothing square. 

Definition 3: contours. Given a spin configuration with q boundary condi- 
tions cr G 6^/ and a rod configuration R E Qa compatible with it, let T be one of 
the maximal connected components of B[a). By construction, the complement of 
r, A\r, consists of one or more connected components: one of these components 
is adjacent to (i.e., it is at a distance 1 from) A"^ and is naturally identified as 
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FIG. 2. The peel of a bad region V 



the exterior of F; it is denoted by Ext F. If F is simply connected this is the only 
connected component of A \ F; if not, i.e., if F has /ir > 1 holes, then there are 
other connected components of A \ F, to be called the interiors of F and denoted 
by IntjF, j' = 1, ■ ■ ■ , h^- The interior of F is then IntF = Ujlnt^F. For what 
follows, it is also convenient to introduce the 1-tile-thick peel of F (see Fig.2): 

Pr = U ^^A': \. (4.4) 

dist'(5,r')=l 

Note that, since distinct /^-disconnected regions are separated by at least one 
smoothing square (i.e., 4 tiles), then also the peels associated to distinct F's are 
mutually D-disconnected. 

The contour J associated to the support F = supp(7) is defined as the collec- 
tion: 

7 = (r, a^, rriext, mm) (4-5) 

where 

• cr^ is the restriction of the spin configuration cr to F; 

• Rj is the restriction of the rod configuration RtoT; 

• ruext is the magnetization of Pp^* :— Ext F n Pr; 

• mm — {''T^im^ • • • ) mm} J with m]„j the magnetization of IntjF fl Pp; if /ir = 
0, then mj^j is the empty set. In the following we shall also denote by 
P^"* := IntF n Pr the internal peel of F. 
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If niext — Q, then we say that 7 is a g-contour. 

Reiiicirk. The set 7 must satisfy a number of constraints. In particular, given 
rriext s-nd ZZlmt; ^7 must be compatible with the conditions that: (i) all the sam- 
pling squares having non-zero intersection with Pp are good (otherwise the con- 
tour would also contain these squares); (ii) each smoothing square contained in 
r has non zero intersection with at least one bad sampling square. Moreover, i?^ 
must be compatible with itself. 

In the following we want to write an expression for Z{A\q) purely in terms of 
contours. Roughly speaking, given a contour configuration contributing to the 
r.h.s. of Eq.(3.6), we first want to freeze the rods inside the supports of the 
contours, next sum over all the rod configurations in the good regions and show 
that the resulting effective theory is a contour theory treatable by the Pirogov- 
Sinai method. The resummation of the configurations within the good regions 
can be performed by standard cluster expansion methods, as explained in the 
following digression. 

Partition function restricted to a good region. Given a set X C A 
consisting of a union of tiles, let = 9 = i- The restricted theory 

of the "uniformly g-magnetized" region X (with open boundary conditions) is 
associated to the partition function: 

Z%X) = Yl ' (4-6) 

which can be easily computed by standard cluster expansion methods, some as- 
pects of which are briefly reviewed here (for extensive reviews, see, e.g., [6] and 
[12, Chapt. 7]). The logarithm of Eq.(4.6) can be expressed in terms of a con- 
vergent series as: 

logZ^X) = Yl z^ViR) = -^1^1(1 + 0{zk)) (4.7) 

where (f^ are the Mayer's coefficients, which admit the following explicit repre- 
sentation. Given the rod configuration R — {ri, . . . ,r„}, consider the graph Q 
with n nodes, labelled by 1, . . . , n, with edges connecting all pairs i,j such that 
i^i n 7^ {Q is sometimes called the connectivity graph of R). Then one has 
(^^(0) = 0, (p^{r) = 1 and, for \R\ > 1: 

^^{R) = ^(_l)number of edges in ^^ g^, 

^- ccg 

where R\ — Hr-^l^)' ^^"^ runs over all the connected subgraphs C of 

G that visit all the n points 1, . . . ,n. In particular, if \R\ > 1, then (p^{R) — 
unless R is connected. 



14 



The sum in the r.h.s. of Eq.(4.7) is exponentially convergent for zk <^ 1; in 
particular, if xq G A, then for a suitable constant C > 

J2 \z\^V{R)\<C<Czkr'^ , (4.9) 

R3xo, \R\>m 

uniformly in A, where R 3 xq means that R contains at least one rod with center 
in xq. Moreover, the sum X^Refj^ • Rg^o ^^^^V^'^i^) analytic in zk, uniformly in A, 
for zk small enough and its limit as A ^ is analytic, too. A useful corollary 
of Eq.(4.9) is the following: if V{R) is the union of the centers of the rods in R, 
supp(/?) is the support of the union of rods r E R (thought of as a subset of A) 
and diam(supp(i?)) is its diameter, then for any finite region X G A: 

V{R)nX^$, R^xo, \R\>m 

diam (supp (il) ) > (i 

<2Cz\X\{Czk)'^^-\ (4.10) 

uniformly in A; here, in the first inequality, we used the fact that in order for 
supp(it!) to have diameter d, the configuration R needs to have at least \d/ (k — l)] 
rods, while in the second inequality we used Eq.(4.9). In a similar fashion, all the 
correlation functions can be computed in terms of convergent series, as long as zk 
is small enough. These results are classical, see [34] or, e.g., [6, 12]. The restricted 
theory is applied to the computation of the sums over the rod configurations in 
the good regions, as described in the following. 

Contour representation of the partition function. Given a contour 7, 
let Zj{lntjr\m.l^^) be the partition function on the j-th interior of F with the 
boundary conditions created by the presence of the "frozen" rods Rj. Moreover, 
if C e P{r, let 

A^(A^) = Ur,ea,{^)\ , C^(Ae) = U ,er'; A, , (4.11) 

dist'(r;,5)<2 

where 

«7(0 {C} U {77 e A' : dist'(?7, = 1, dist;(?7, T') = 2, r)j^_,) = , (4.12) 

with disti(-, •) the rescaled ("coarse") Li distance on A' and — 1, — 2. 
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FIG. 3. The two sets A^(A^) and C-y(A^) in the case that q = +. 



Finally, given A C Pr, let /a and qa be the following characteristic functions: 

1 if i? has at least one rod belonging to ^7 (A) 
/a(-R) = •{ and one belonging to (^^(A) , (4-13) 

otherwise , 

1 if i? n -R-y 7^ 0, -R has at least one rod belonging to A^{A) 
gA(R) — { and has at least one rod belonging C-y(A) , 

otherwise . 

(4.14) 

Pictorially speaking, /a is the characteristic function of the event "R crosses the 
boundary of F at A", while qa is the characteristic function of the event "i? 
intersects R^ across A". Note that, by construction, given two distinct tiles, 
Ai C Pn and A2 C Pp^ such that Ai n A2 = 0, then ^^^(Ai) n ^^^(As) = 0, 
even in the case that Fi = F2. 

In terms of these definitions, the following contours' representation for Z{A\q) 
is valid. 

Lemma 1 The conditioned partition function Z{A\q), q — ±1, can be written as 

Z{A\q) = Z'^iA) Yl [U^,il)] ^'"^^'^ , (4-15) 
deC{A,q) -red 

where: 

• C{A,q) is the set of all the well D- disconnected q-contour configurations 
in A (here we say that {71, . . . ,7^} is well D-disconnected if the supports 
Fi,...,F„ are separated among each other and from A^ by at least one 
smoothing square); 
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Cqil) is the activity of ^: 

Ul) = C;(7) cxp { - 5^ ^^{R)z\^\ F^{R)} , (4.16) 



where 

and Fa = /a A C Pf"* w/iz/e Fa = /a + ^a(1 - /a) «/ A C F^^*. 

• is the interaction between the contours in d: 

* 

W(5)= J] ^W'''E(-1)"^' E ^a,(F)---Fa„(F), (4.18) 

RgCII ri>2 Ai<-<A„ 

where the * on the sum indicates the constraint that Ai,...,A„ are all 
contained in the peel of some contour of d and their centers ^i, . . .^n all 
belong to the same row (if q = +) or column (if q = —) of A' , namely 
Ci,j{-q) — ■ ■ ■ — ^n,j{-q)- Morcovcr, by writing Ai < • • • < A„, we mean that 
< ■■• < Cn,j(g)- Finally, Fa = /a ^/ A is contained in the internal 
peel of some contour in d or Fa = /a + 5'a(1 — /a) ^/ A is contained in the 
external peel of some contour in d. 

Remarks. 

1. The contour configurations {71,..., 7„} G C(A,g) consist of n-ples of 
well D-disconnected g-contours, which means that the geometric supports 
Fi, . . . ,r„ are separated among each other and from by at least one 
smoothing square. Note, however, that their external and internal mag- 
netizations are not necessarily compatible among each other: for instance, 
Fi may have one hole surrounding F2, and the internal magnetization of 
Fi may be different from the external magnetization of F2 (which is q). 
It is actually an important point of the representation Eq.(4.15) that we 
can forget about the compatibility conditions among the internal and ex- 
ternal magnetizations of different contours. There exist different (and even 
more straightforward) contour representation of Z{K\q) where the inter- 
nal and external contours' magnetizations satisfy natural but non-trivial 
constraints (e.g., in the example above, the natural constraint is that the 
internal magnetization of Fi is the same as the external magnetization 
of F2). However, the magnetization constraints are not suitable to ap- 
ply cluster expansion methods to the resulting contour theory. Therefore, 
it is convenient to eliminate such constraints, at the price of adding the 
extra factors Z^(IntjF|m]^J/Z(IntjF|g) in the definition of the contours' 
activities, see Eq.(4.17). 
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2. The interest of the representation Eq.(4.15) is that the contour activities 
and the multi-contour interaction satisfy suitable bounds, allowing us to 
study the r.h.s. of Eq.(4.15) by cluster expansion methods. In particular, 

sup;^ Eij^e^rK) 1^9(7)1 < exp{-(const.)2;A:2|r|}, where the * on the sup 
reminds the constraint that all the smoothing squares in F must have a non- 
zero intersection with at least one bad sampling square. Moreover, W{d) 
is a quasi-one-dimensional potential, exponentially decaying to zero in the 
mutual distance between the supports of the contours in d. The proofs of 
these claims will be postponed to the next sections. 

Proof of Lemma 1. Given a e Q\, a spin configuration with q boundary 
conditions consider the corresponding set of contours {71, . . . , 7„}. Some of them 
are external, in the sense that they are not surrounded by any other contour in 
{71,..., 7„}. By construction, these external contours are all g-contours. We 
denote by Cexti^,(l) the set of external g-contour configurations. Given d G 
Cext{^,(l), there is a common external region to all the contours in d, which we 
denote by Ext (9). Besides this, there are several internal regions within each 
contour 7 G 9. For each external contour 7 G (9, we freeze the corresponding 
rod configuration and sum over the rod configurations inside all the internal 
regions Int^F, j = 1, . . . , hr- In this way, for each such interior, we reconstruct 
the partition function Z^(IntjF|m]^j). On the other hand, by construction all 
rods inside Ext(9) are either horizontal or vertical, according to the value of q. 
Therefore, if we sum over all the allowed rod configurations inside this region we 
get the restricted partition function Z|(Ext(9)), where the subscript d reminds 
the fact that the rods Rq — U^^g create an excluded volume for the rods in 
Ext (9). Using these definitions, we can rewrite 

hr 

Z{A\q) = ^a(Ext(a)) n [piR,) n Z,{Int,r\mij] . (4.19) 

Note that here we used the fact that the exterior and the interior(s) of d do not 
interact directly (i.e., they only interact through R^). Using the definition of 
^^(7), Eq.(4.17), we can rewrite Z{A\q) as 



E n[c?wn^]e-'""--. (-0) 



where 



-ws-O) _ ^a(Ext(a)) U,ea Ujii ^^Int^F)] 

~ Zi{A) ■ ^ ' 

The factors ^Jq^intT) same form as the l.h.s. of Eq.(4.20) itself, with A 

replaced by Int^T: therefore, the equation can be iterated until the interior of all 
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the contours is so small that it cannot contain other contours. The result of the 
iteration is 

Z{K\q) 



Zi{K) 



E [nC(7)]e-^°(^), (4.22) 



where 



Zi(A) ' ^ ' ' 

A(9) = A\ U^ffzQ r is the complement of the contours' supports and Zg(A((9)) is 
the restricted partition function with magnetization q in the volume A{d) and in 
the presence of the hard rod constraint generated by the frozen rods R.y in the 
region U^gg U^cp^^* ^7(^)- 

We now use Eq.(4.7) and the analogous expression for Zg(A(9)), i.e., 



ext 

R n Rg=tD 

logZ|(A(c?))= Yl ^^V{R), (4.24) 



ext 

where Rn Rg = means that R does not intersect Rg from the outside, namely: 



RnRg^$ ^ n n (l-^A(it:)) = l, (4.25) 

where g^A was defined in Eq.(4.14). Then we can rewrite: 

,-Wn(a)^ ^^(^(^))n,ea^^(n Z^(A(g)) 

Z^iA) Zi{A{d)) ' ^ ■ ^ 

= exp { - 5^ z\\^{R)} ■ exp { - 5^ z\\^{R)} , 



g 

where i?-^2 means that R must contain two rods ri,r2 belonging, respectively, 
to two distinct elements of the partition V{d) of A induced by the contours in 
d; i.e., either ri,r2 G i? belong, respectively, to two disconnected components of 
A{d), or they belong to two different contours' supports, or ri belongs to one 
contour's support and r2 to one of the components of A{d). Using the definitions 
of the characteristic functions /a and qa defined in Eqs.(4.13)-(4.14), the two 
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exponential in the r.h.s. of Eq.(4.26) can be written as 

X; z\^V(R) - E -""^^(^) [l - n 11(1 - MR))] , (4.27) 

ReniiA) Real jedAcPr 

E z^ViR) = E ^^ViR) [l - n n (1 - 9a{R))] ■ (4.28) 

i?en« Real -yed acp=^' 

•[n 11(1 -/A(i^))" . 

Using the representations Eqs.(4.26), (4.27), (4.28) into Eq.(4.22), we find 

E [nC(7)]exp{-E-'VW- (4.29) 

^ ^ deC{A,g) -red Renl 

•[i-(n n (i-^A(i?)))-(n n(i-/A(i?)))]}- 

Note that the expression in square brackets in the second fine can be conveniently 
rewritten as 

i-n( n (i-5A(i?))(i-/A(i?)))-( n (i-/A(i?)))^ 

7G9 ACP=^* ACPf"* 

^1-n n(i-^A), (4.30) 

769 ACPr 

where Fa was defined in the statement of Lemma 1. Plugging Eq.(4.30) into 
Eq.(4.29) gives 

= E [nC(7)] exp{ - E ^""V'^W ■ (4-31) 

^ aeC(A,g) 7ea Renl 

•[E^A(i?) + E(-ir^' E ^A,(i?)---FA„(i?)] , 

ACPa n>2 {Ai,..,A,J 

where the sum a„} I'ui^s over collections of distinct tiles Aj C Pq. with 

:= U^ggPr- Finally, using the fact that ^'^{R) forces R to be connected and, 
therefore, to live on a single row or column, depending on whether g is + or — , we 
find that the only non- vanishing contributions in the latter sum come from n-ples 
of tiles all living on the same row or column. This proves the desired result. ■ 



5. REORGANIZING THE CONTOUR EXPANSION. 

Standard cluster expansion methods are more easily implemented in the case 
of two-body interactions. Our contour interaction Eq.(4.18) is many-body but it 
can be reduced to the two-body case by a slight reorganization of the expansion. 
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Lemma 2 The contour representation (4-15) for the conditioned partition func- 
tion Z{A\q), q — ±1, can be reorganized as follows 

^^ = i+E E (^i)---i^fn^^)0({^i,^2...,x„}), (5.1) 

^ ^ m>l{Xi,...,Xm} 

where: 

• each polymer Xi is a D-connected union of tiles in A; 

• (f) implements the hard core interaction, i.e., 



(f){Xi,Xj) 



i<j 

1 if X, D- disconnected from Xj 
otherwise. 



is a suitable function of X, called the polymer's activity, which is 
defined by Eq.(5.3) below. 

Remark. The definition Eq.(5.2) of the polymer interaction is the analogue 
of Eq.(3.2) with the rods replaced by polymers and the notion of intersection 
replaced by D-connectedness. 

Definition of the polymer's activity. Given d = {71, . . . 7„} G C(A, q) and 
Xq = U"^^rj, let Y = {Ai,. . .Am} be a collection of m > 2 distinct tiles, all 
contained in the peel of Xq, i.e., Aj G U^^^Prj, and all belonging to the same 
row (if g = +) or column (if q = —)■ Since the tiles are all on the same row 
(column), we can order them from left to right (bottom to top), Ai < A2 < 
• • • < Am- We denote by T^^ the set of all such collections. Moreover, for each 
Y — {Ai, . . . A^} G with Ai < • • • < A^^, we define Y to be the union of all 
the tiles between Ai and A^. With these definitions, the activity of the polymer 
X is given by 

E E [nc.(7)][n(^^^"'^-i)]' (5-3) 

n>l,p>0 aeC(A,g): \d\=n jed i=l 



XsU{UjYj}=X 



where Cgil) was introduced in Eq.(4.16) and, if y = {Ai, . . . , A^} G T^^, 

HY) ■= (-ir E z^'^VrnFAM-'-FAAR) ■ (5.4) 
Reni 



21 




FIG. 4. An example of polymer X and of a possible way of realizing it as a union of 
three contours' supports Fi, F3 and of two sets Yi, Y2. 



An example of a polymer X with non- vanishing activity and of a possible way of 
realizing it as a union of sets Fj and Yj is given in Fig. 4. 

Proof of Lemma 2. Using the definition of J^{Y), we can rewrite Eqs.(4.15)- 
(4.18) as: 



Zi{A) 



= E [new 

aeC(A,g) 7ea 



n 



MY) 



(5.5) 



Let us now add and subtract 1 to each of the factors e-^^^\ In this way we turn 
each factor into a binomial 1 + {e-^^^'> — 1). If F G '^Xg^ associate the quantity 
(^Q-^(Y) — 1) with the region Y; similarly, we associate the activity (^(7) with the 
region F. In this way, every factor of the form IliLi C(7j) 11^=1 (^"^'■'^^■' ~ 1) 
geometrically associated with the region X = {U"^;^Fj} U {U^^{Yj}. We develop 
the binomials 1 + (e-^*^^) — 1) and collect together the contribution corresponding 
to the maximally D-connected regions, obtained as unions of Fj's and Yj^s. The 
result is equation Eq.(5.1). ■ 



6. CONVERGENCE OF THE CONTOURS' EXPANSION 

In this and in the next section we prove the convergence of the cluster ex- 
pansion for the logarithm of the partition function with q boundary conditions, 
starting from Eq.(4.15). The proof will be split in two main steps: first, in this 
section, we prove convergence under the assumption that the activities Cgil) sat- 
isfy suitable decay bounds in the size of |F|. Then, in the next section, we prove 
the validity of such a decay bound via an induction in the size of |F|. Prom now 
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on, C, C", . . . and c, c', . . . indicate universal positive constants (to be thought of 
as "big" and "small", respectively), whose specific values may change from line 
to line. 

Lemma 3 Suppose that, for zk and {zk'^y^ small enough, 

sup* Ul)\<e-'^'''^^\ (6.1) 



where the * on the sup reminds the constraint that all the smoothing squares in 
r must have a non-zero intersection with at least one had sampling square, and 
Cq — b* 10~^. Then the logarithm of the partition function admits a convergent 
cluster expansion 

logZ(A|g) = z\%\R) + E [ n K\^\^f{X) , (6-2) 

where X — {X\, . . . , X„} is a polymers' configuration (possibly, some of the Xi 's 
may coincide), each polymer X being a D-connected subset of A consisting of a 
union of tiles. 

Remarks. 

1. The constant Cq — 5 * 10~^ is a possible explicit constant for which the 
result of the lemma holds (certainly, it is not the sharp one). Its specific 
value is motivated by Lemma 4 and by its proof, see next section. 

2. The function in Eq.(6.2) is the Mayer's coefficient of A", defined as 
in Eq.(4.8), with R replaced by A", by X^, and the notion "rj fl rj ^ 0" 
replaced by "Xj is D-connected to Xj" . 

Proof. By Lemma 2, Eqs.(4.15)-(4.18) can be equivalently rewritten as 

§^ = i + E E (^™)0({^i,---^m}). (6.3) 

^ ^ m>l {Xi,...,X^} 

It is well-known [6, 12, 34], that if the activities Kq^\x) are sufficiently small and 
decay fast enough with the size of X, then one can apply standard cluster expan- 
sion methods (analogous to those sketched above, after Eq.(4.6)) for computing 
the logarithm of Eq.(5.1) and put it in the form of the exponentially convergent 
sum. More specifically, a sufficient condition for the application of the standard 
cluster expansion is, see e.g. [12, Proposition 7.1.1], 

|Kf)(X)|<C£|f'le-'^°^'(^'), (6.4) 
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for some kq > and Sq small enough (here 5'{X') is the rescaled tree length of 
the coarse set X' C A', i.e., it is the number of nearest neighbor edges of the 
smallest tree on A' that covers X'). In the following, we will prove that under 
the assumption of the Lemma, the polymers' activities Kg^^ (X) satisfy 

\Kjj^\X)\<e,£\'''\-\ £i:=e->^% := (zk)^^ , £ := max{£i, £2} , 

(6.5) 

where cq is the same constant as in Eq. (6.1). Using the fact that X is D-connected, 
we see that Eq.(6.5) implies Eq.(6.4) with £0 = e'^o = e^/^. Therefore, by [12, 
Proposition 7.1.1], we get 

i°sTi = E[]lK'"m]nx). (6.6) 

^ ' xcA xex 

Combining this equation with Eq.(4.7) gives Eq.(6.2). 

The rest of this section is devoted to the proof of Eq.(6.5). The polymer's 
activity Eq.(5.3) can be rewritten as 

<'m= EE E [iluii)]- (6.7) 

Xo,XiCX: n>l a={7i,...,7„}eC(A,g): j=l 

XouXi=x u^^iri=Xo 



n(e^«) - 1) 



.1=1 



■EE E 

UiSupp{Yi)=Q 

Note that since all the tiles in Q belong to the peel of some contour then Q nXo = 
0. On the other hand, the sets Xq and X^ may very well overlap fl 7^ 
in general (see Fig.4 for an example). Moreover, once Xq is fixed, the supports 
of the contours are automatically fixed too, since they must be the D-connected 
components of Xq. Then we can rewrite the sum as 

m^)- E E [nc.(7.)]- (6.8) 

07^XoCX {7i,...,7„}eC(A,g): J=l 
supp(7i)=ri(Xo) 

• E EE E [n(^^™-i) 

XiCX: QCXi p>0 {Yi,...,Yj,}: YiST^ i=l 
^"''^1=^ UiSupp(yi)=Q ° 

UiYi=Xi 

where in the second sum Ti^Xo) are the maximally D-connected components of 
Xq, which must be well D-disconnected (otherwise the corresponding contribution 
to the activity is zero). 
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Now, note that T{Y) given in Eq.(5.4) is at least of order n (with n > 2) 
in z, by the very definition of the characteristic function Fa. In fact, F^{R) 
is either equal to /a or to /a + 5'a(1 — /a); therefore, using the definitions of 
/a and (7A) Eqs.(4.13)-(4.14), we see that F^{K) is different from zero only if 
R contains a rod belonging to A^(a)(A). Now recall that, as already observed 
after Eq.(4.14), distinct tiles Ai 7^ A2 correspond to distinct sets A^(Ai)(Ai) and 
A^(A2)(A2), such that A^(Ai)(Ai) fl A^(A2)(A2) = (here 7(Aj) is the contour 
whose peel Aj belongs to, Aj G Pr^ since the peels of different contours are 
disconnected, the contour 7(Ai) is unique). Therefore, the r.h.s. of Eq.(5.4) is 
non zero only if R contains at least n distinct rods. Using Eq.(4.9), we find that, 
if y = {Aj,, . . . , A^„} with A^, < • • • < A^„, 

|7-(F)| < J2 \4V{R)\ < 2Czf{Czkr^^'^-'''^ 

Renl: \R\>2 
y(R)nAi^0, 

diam(supp(il)) >diam(y ) 

< C'zeizky^'^'^'^'"^ , (6.9) 

where diam'(y) = — is the rescaled diameter of the set UAeyA, and a 
can be chosen to be a = 1/4. Using this bound and the fact that |e^ — 1| < |j;|el^l, 
we find: 

J2 n|e^«)-l|< (6.10) 

{Yi,...,yp}: YiST^^ i=l 
UiSupp{Yi)=Q 
UiYi=Xi 

< fl {^^'^^^(^^)"-^'^"^'«)e[^'^'=^(^'=)"-^'^"^''1 1 

{Yl,...,Yp}: Yi&T"^^ i=l ^ 
UiSupp(Yi)=Q 

UiFi=Xi 

Now, since the choice of Y only depends on the union of the contours' supports 
Xq, in Eq.(6.8) we can start with performing the sums over the contours' spin 
attributions, rod configurations and internal colors. Using the bound Eq.(6.1) on 
the contours activities, we get 

n n 

E [nc^M^n[ E (6-11) 

{7l,...,7n}eC(A,g): j=l j=l 7i- 

supp(7i)=ri(Xo) supp(7j)=rj(Xo) 

<Jj6l^^le"'^°'^'='l^^l = Q\xi>\^-co-zk^\xi,\ 
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where the factor 6 1 '"'I bounds the sums over cr^ and m^^^ at F fixed. Putting these 
results together into Eq.(6.8), we find 

9j^XoCX XiQX: QCXi p>0 {Yi,..,yp}: FieT^^ 

XoUXi=X , 

UiYi=Xi 

■ n |c'.)t^(.)t)-^^-'«)e[^'^'=^(^^)"""'"''^1 1 , (6.12) 
j=i ^ J 

Now, note that: (i) ^J^^ diam'(Fj) > |X(| - 1 > \X[\/2; (ii) > \Q'\, 

because every tile in Q = Uisupp(l^) belongs to the peel of Xq; 

p _ 

(iii) E(^'^)""^'^"'^'''^ ^ E E^^^)"'""'' - c"\Q'\{zk)'' 

Plugging these estimates into Eq.(6.12) we find 

|i^f)WI< E (6e-'^"'')l''°l E (^^)^"'^'- (6-13) 

$^XoCX XiCX: 
XoUXi=X 



Qcxi p>o Yi&r" j=i 



•^0 

LliSupp{Yi)=Q 
UiYi=Xi 



which can be further bounded by: 



\K^^\X)\< E e-*^^'l^ol J2 (6.14) 

$^XoCX XiQX: 
XoUXi=X 

Qcxi p>o ^' AnQ^0 

|A'|>2 

where in the last sum A is a generic subset of A consisting of a union of tiles, and 
S'{A') is its rescaled tree length. The expression in square brackets in the second 
line is bounded above by C"zk'^\Q'\{zk)^ , so that 



0^XoCX XiCX: QCXi 

XoUXi=X 



zk^\Q'\{^-C"{zk)^) 



< J2 e-^^^'l^ol J2 (^^)^"'^' E e-^^^'l^'l . (6.15) 



$^XoCX XiQX: QCXi 

XoUXi=X 
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The last sum can be rewritten as X^qcxi ^ "^'^^''^'l = (1 + e 4 so that, 

defining Si :— e~^^^ , £2 '■— (zk)^ and £ := max{£i,£2}: 

X(,UXi=X 

< Y: 4^^' E (2^2)1^^! = (ei + 2£V~2 + 2£2)l^'l-(2£2)l^'l 

H^XoQX XiCX: 
XoUXi=X 

< \X'\ei{l + 2^2) (£1 + £i2£2 + 2£2)'''''"' < (6.16) 

where e — max{£i,£2}- Setting ei = y/ii, £2 = and recalling that a = j, we 
obtain the desired estimate on Kq^\x). This concludes the proof of the lemma. 
■ 

Remcirk. The dependence of the activities Kq^\x) on A is inherited from 
the constraint that X must be separated from A'^ by at least one smoothing 
square, and by the fact that the quantities ^(t) and -FiY) themselves are A- 
dependent, simply because their definitions involve sums over rods collections in 
fl\. However, this dependence is very weak: in fact, if Kq{X) is the infinite 
volume hmit of Kq^\x), we have: 

\K(^\X) - KqiX)\ < (^,el^'l-i)'V"i^^*'(^'<) , (6.17) 

for some c' > 0. The proof of Eq.(6.17) proceeds along the same lines used to 
prove Eq.(6.5) and, therefore, we will not belabor the details of this computation. 



7. THE ACTIVITY OF THE CONTOURS 



In this section we prove the assumption Eq.(6.1) used in the proof of Lemma 
3. Let us first remind, for the reader's convenience, the definition of Cq(7)- 

C,(7) = C;(7) exp { - ^ ^^{R)z\^\ ^A(i?)} , (7.1) 

where 

^1^^' zi{r) j}^ z(int,r|g) ■ ^ ' 

By using the same considerations used to get the bound Eq.(6.9), we see that the 
expression in braces in the r.h.s. of Eq.(7.1) is equal to a contribution of order one 
in z plus a rest, which is bounded in absolute value by Cz/c^ |r'| (zk)'^. On the other 
hand, the contribution of order one in z is equal to —z X^^g^?. X^ac/v -^a(-R), 
which is negative, simply because Fa > 0. Therefore, 

lC.(7)l<lC?(7)|e^^'^^l^'l(^'=)", (7.3) 
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which makes apparent that, in order to prove Eq.(6.1), we need to prove an 
analogous bound for Cqil)- By definition, Z-y{X\m) < Z[X\m), so that 

icwi.ic:win-{i.M_^}. c:w:^|^. (7.4) 

The estimate that we need on the quantities Cg(7) ^-iid ^zQ^ii T\q) '^^ summarized 
in the following two lemmas. 

Lemma 4 Let zk and {zk'^)~^ be small enough. Then 

sup* E IC;(7)l<e-^-^'=^l^'l, (7.5) 



^ R^eOr(o"7) 



where the * on the sup reminds the constraint that all the smoothing squares in 
r must have a non-zero intersection with at least one bad sampling square, and 
Co = 5 * lO"! 

RemEirk. The specific choice of cq in the lemma comes from Eq.(7.10) below. It 
is related to the size of the smoothing squares, to the number of zero spins and 
to the number of pairs of neighboring spins with opposite sign that can appear 
in a contour (as explained below, it comes from the remark that every smoothing 
square - which contains 64 tiles - in a contour must intersect at least one bad 
samphng square - of size £^ > 

Lemma 5 Let zk and {zk'^)~^ be small enough. Then there exist two positive 
constants C, Ci > such that, for any simply connected region X dl? consisting 
of a union of smoothing squares, 

-\P'^\{Cze{zk)+e-i) < ^(^1 + ) < \P'^\{Czk\zk)+e-i) (n 

-Z{X\-)~ ' ^ ' 

where e was defined in Eq.(6.5) and Px is the 1-tile-thick peel of X. 

These two estimates combined with Eq.(7.3) give 



2|r'l r>'^^.2|^'l/,^^" TT \P{^^ ,p\(C zk^ (zk)+e''i) 



sup* K.(7)l< e-^'^'^^'\^'\e^'^^'\^'\^^^^" JJe 



(T'y 



< ^\T'\{Czk^(zk)+e''i)^C'zk''\r'\{zkr^-2co zk''\r'\ 



^ ^-zk''\r'\{2co-c'{zkr-c(zk)-^) ^ ^-cozk'^\r'\ 

under the only assumptions that zk and {zk"^)'^ are small enough. Therefore, 
these two lemmas imply the convergence of the cluster expansion Eq.(6.2), which 
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completes the computation of the partition function of our hard rod system with 
q boundary conditions. A computation of the correlation functions based on a 
similar expansion will be discussed in the next section. The rest of this section 
is devoted to the proofs of Lemma 4 and 5. 



Proof of Lemma 4. Let be a spin configuration compatible with the 
fact that 7 is a contour. In particular, let us recall that every smoothing square 
contained in F has a non zero intersection with at least one bad sampling square; 
moreover, by its very definition, each such bad square must contain either one tile 
with magnetization equal to 0, or one pair of neighboring tiles with magnetizations 
+ and — , respectively. Therefore, given cr-j,, it is possible to exhibit a partition 

of r such that: (i) all the elements of the partition consist either of a single 
tile or of a pair of neighboring tiles with opposite magnetizations + and — (we 
shall call such pairs "domino tiles" ) ; (ii) if is the number of single tiles in V 
with magnetization equal to and Ha is the number of domino tiles in then 

-\- Md > I r' 1/64. The factor 64 comes from the consideration that in F, by 
definition, we have at least one bad square every four smoothing squares, and by 
the fact that four smoothing squares contain 64 tiles. 

By the definition of (p{R^), we have: (^{R-y) < Wp^-p '~p{Rp) ■ Moreover, using 
the standard cluster expansion described after Eq.(4.6), we find that Z^(F) > 
ripe-p Z^(P)e~'"^^^'^^^^l^'l. By combining these two bounds we get 

E ic:(7)i<e°'"<""^''n|EiMl' 

where the sum over Rp runs over rods configurations in J7p(U^gp/(T^). Now, if P 
is a single tile with magnetization either + or — , then YliRp ^i^pj = 1- Moreover, 
if P is a single tile with magnetization equal to 0, then J^^^ '^z^p) ~ ~z^p) ~ 

_^~zP{\+0{zk)) _ 

Finally, let us consider the case that P is a domino tile. We assume without 
loss of generality that P = {A^j, A^j}, with ^2 — ^1 = (^, 0), and cr^^ = "Cr^ — +• 
Since the rods interact via a hard core, ip{R^-^, R^^) is different from zero only if 
at least one of the two rod configurations R^^ and R^^ is untypical: here we say 
that P^-^ is untypical if it does not contain any rod in the right half of A^^ and, 
similarly, that R^^ is untypical if it does not contain any rod in the left half of 
A^j. Therefore, 

r fiR^i) , ^{R^2) 1 

R(^ untypical R^^ untypical 
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where we used that X^Ren^ fi^) — -^^(A). Eq.(7.9) can be rewritten and esti- 
mated (defining A|'^ to be the left half of A^J as 

'fjRp) < n„Czk^(zk) ^^^HJ < 2e^'^*^'(^'=)e-^^'/^ (7 10) 

^ Z<i{P) - Z+(A^J - ■ \ - ) 

Plugging the bounds on i^^to Eq.(7.8) gives: 

< g-^£2(l-C'^fe)|r'|/128 ^ 

where in the last line we used the bound A/q + Md > |r'|/64. Using I > k/2 we 
obtain Eq.(7.5) so the proof of the lemma is complete. ■ 

Proof of Lemma 5. We proceed by induction on the size of X. If X is so 
small that it cannot contain contours D-disconnected from X'^, then 



Z{X\+) _ Z+{X) 

z{x\-) - z-{x) 



exp{ ^ - . (7.12) 



Let V{R) be the union of the centers of the rods in R and let R G fi^. Since the 
orientation of all rods in R is fixed, V{R) identifies uniquely the rod configuration. 
Then 

sex ReQ^2 xeXi?en«2\f2x 



The first sum in the second line is equal to 

)^(i?)zl- 

\vm 



where 

\v{R)\ 



V(R)3x 



is an analytic function of z, of the form s{z) — z{l + 0{zk)), independent of q 
and X. The second sum in the second line of Eq.(7.13) involves rod configurations 
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containing at least one rod belonging to X and one belonging to X'^. Therefore, 
it is of order at least 2 in z and scales like the boundary of X: 

E E ^SMr|sc..~*'^(-*)i^'ii. (7.16) 

V{R)3x 

for a suitable constant Ci > 0, independent of q. Plugging Eqs.(7.13)-(7.16) into 
Eq.(7.12) gives: 

Z{X\-)~''''^\ ^ ^ \V{R)\ +^ ^ \V{R)\ /'^^--^^^ 

y(R)9a; V{R)3x 

which is bounded from above and below by e'^Cizk'^{zk)\p'x\ ^-2Cizk^(zk)\p^\ ^ 
respectively. Setting C > 2Ci, this proves the inductive hypothesis Eq.(7.6) at 
the first step, i.e., for regions X small enough. 

Let us now assume the validity of Eq.(7.6) for all the regions of size strictly 
smaller than Aq, and let us prove it for Aq. As explained in Section 6, Z{Ao\q) ad- 
mits the cluster expansion Eq.(6.2) involving polymers X that are D-disconncctcd 
from Aq, whose activities arc defined in Eq.(5.3). In particular, the cluster ex- 
pansion is convergent provided that Cij(7) is bounded as in Eq.(6.1). Now, note 
that the interiors of the contours 7^ involved in the cluster expansion for Z{Xo\q) 
via Eqs.(6.2) and (5.3) have all sizes strictly smaller than Aq. Therefore, us- 
ing the inductive hypothesis, the product max |l, ^^(int^rjq) } Eq.(7.4) can be 
bounded from above by e^'^'^^'-^^'^'^^^'^^'^^"'' that, if combined with Eqs.(7.3), (7.5), 
implies Eq.(6.1) for all the the contours 7^ involved in the cluster expansion for 
Z(Xo|g). We can then write: 

= §^-P { E . (7.18) 

where Ki^°\x) = Uxex and Ki^°\x) admits the bound Eq.(6.5). 

The first factor in the r.h.s. of Eq.(7.18) is rewritten as in Eq.(7.17) and is 
bounded from above and below by e'^^^^'^^^^'^^^^^o^ and e~'^^'^''^'^^''^^^^o \ respectively, 
exactly in the same way as Eq.(7.18) itself. 

The second factor in the r.h.s. of Eq.(7.18) can be bounded as follows. We 
rewrite 



exp { J2 - K^^°\x)]cl>'^{X)} = (7.19) 

XCAo 

= exp{ J] 5X,(A')0^(A')} -expl J] g[i^f°)(A') -X,(A')]0^(A')} , 



XCAo XCAo 
q=± q=± 
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where Kg{X) = U^^^ Kq{X). Now 



n n n m—l n 

j=l j=l m=l j=l j=m+l 

(7.20) 

then using Eq.(6.17) and (6.5) , we have 

n 

|Xf °)(A') - K,{X)\ < ^(^)l^;.l£c'dist'(x;„Ai,j j-j ^\x'.\ 

m=l 

< ^c'dist'(^,Ai,J [J £^ , (7.21) 

Therefore, the second factor in the second hne of Eq.(7.19) can be bounded from 
above and below by e ' '^o' and e ' ^o', respectively for a suitable constant C2- 
We are left with the first factor in the second hne of Eq.(7.19), which involves 
the partition sum 

Y, K,(XW(X) ""''Ivf''''' . (7.22) 



XCAo 



where [A"'! is number of tiles in UxexX. Eq.(7.22) can be further rewritten as 



XCAo ^eA'o XDA^ I I 

where 

s...j:'SMm^, (,,4) 



XDA^ 



XCl? 

is independent of q and ^. The second term in the r.h.s. of Eq.(7.23) is bounded 
in absolute value from above by I-PaoI^'^^ ^ suitable C3 > 0; therefore, 

exp { iK,(XW(X)\ = exp { ^ ^ i"^^^^^] < e^''"-"" 

XCAo ^eA'o XDA^ I I 

9=± q=± XnA^^$ 

(7.25) 

and is bounded from below by e '^o . Choosing c\ such that e'^^ > e^^ + 2e^^ 
this completes the inductive proof of Eq.(7.6). ■ 



8. EXISTENCE OF NEMATIC ORDER 



In this section we prove Theorem 1. We start by proving Eq.(3.9). The 
probability that the tile centered at ^0 has magnetization —q in the presence of 
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boundary conditions q can be written as 

d 



{x^o%-i^^ogZM<l) , (8.1) 

OZq zo=1 

where Zzo{A\q) is defined in a way completely analogous to Eqs.(3.4)-(3.6), with 
the only difference that the activity C(0 Eq.(3.4) is replaced by C(0) where 
m = C(0 if e 7^ 6, while 

z\^(\ if = g 
C(Co) = <; Zoz\^^\ iia^ = -q (8.2) 
-1 if (7^ = 0. 

The change of C(0 i^^^o C({) induces a corresponding change of Cg(7) and ^^^''(X) 
into C(j(7) ^iid Kq^\x), respectively. The activity Kq^\x) admits the same 
bound Eq.(6.5) (possibly with a slightly different constant c"), uniformly in Zq 
for Zq close to 1, and it depends explicitly on Zq only if X 3 A^(,. In such a case, 
the derivative of Kg^\x) with respect to zq is bounded by ^/ eis^^'^~^, uniformly 
in Zq for Zq close to 1. 

The logarithm of the modified partition function admits a convergent cluster 
expansion analogous to Eq.(6.2): 

^ ' XCA 

SO that 

(X^o')l = E ^-^f W ■ (8.4) 

XCA 

The sum in the r.h.s. of Eq.(8.4) is exponentially convergent for e small enough, 
and it only involves polymer configurations containing A^^, simply because 
Kq^\x) is independent of Zq whenever A^^ n X = 0. Therefore, 

XCA ^^^(o 

< (const.) (£i) 2 , (8.5) 

which proves Eq.(3.9). 

In order to compute the density-density correlation functions we proceed in a 
similar fashion. We replace the activity 2; of a rod r centered at x by z^ and we 
define Zz{A\q) to be the modified partition function with boundary conditions q 
and variable rod activities z = {zx}xga- Correspondingly, we rewrite: 

{nx)l = zd,JogZ^{A\q) 

z=z 

- Ml{ny)l = z^d,J,^ logZ^(A|g) , (8.6) 
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where — z means that = Vx € A; the higher order density correlation 
functions have a similar representation. Once again, logZz(A|g) admits a cluster 
expansion completely analogous to logZ(A|g): 

log^z(A|g) = [n^-w]^^(^) + E ^ifzH^)0^W| , (8.7) 

where x{r) is the center of r. Moreover, Kq^^{X), together with its derivatives 
with respect to and/or Zy, admit the same bound Eq.(6.5), possibly with a 
different constant c"; the derivative of Kq^^{X.) with respect to Zx and/or Zy is 
different from zero only \i X 3 x and/or X 3 y. Therefore, 

{nxny)l - {nx)l{ny)l = 

= Y z\^^R{r{x))R{r{y))^^{R) + Ydl.K^q^{X)<P^iX) 

where R{r) is the multiplicity of r in R. The sums in the first line involve 
connected rod or polymer configurations containing at least one rod centered 
at X] similarly, the sums in the second line involve connected rod or polymer 
configurations containing at least one rod centered at x and one rod centered at 
y. All the sums are exponentially convergent and their evaluation finally leads to 
the finite volume analogues of Eqs.(3.10)-(3.11). The infinite volume counterparts 
are obtained simply by replacing all the finite volume activities with their infinite 
volume counterparts and by dropping the constraints that the polymers should be 
contained in A. The infinite volume limit is reached exponentially fast and all the 
observables share the same invariance properties as the infinite volume activities 
themselves. In particular, the infinite volume Gibbs measures (•)^ are translation 
invariant, and the averages {x^^Y ^'^d (W^n^^Y are all independent of q. We 
will not belabor the proofs of these claims, since they are all straightforward 
consequences of the cluster expansion described in the previous sections, in the 
same sense as the representations for (x^'^)^, {uxY and {nxUyY and the proof of 
their convergence, discussed in this section, are a consequence of the bounds of 
sections 4 and 6. This concludes the proof of the main theorem. ■ 
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